Ab initio investigation of the groundstate, magnetic, electronic, and optical properties 

of polyyne and cumulene prototypes 
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We have investigated polyyne and cumulene prototypes based on the density-functional theory. 
Our independent-particle spectra show that the various carbynes can be distinguished by optical 
properties comparing the low-energy spectral structure as well as using very general considerations. 
The latter conclusion is supported by results based on the random-phase approximation including 
local-field effects. 

PACS numbers: 31.15.A- 31.15.es 78.67.-n 31.15.ap 



I. INTRODUCTION 

Carbon-based materials possess many potential ability 
for future technologies.*^ Carbon atoms (C) are known 
to have three bonding states resulting from hybridization 
of the atomic orbitals: sp 3 (diamond) , sp 2 (graphite and 
fullerenes), and sp. While the importance of sp 2 struc- 
ture is by now confirmed, a lot of interest is growing 
around sp carbon allotropes, named carbynes. 

The formation of sp carbon chains is expected in 
the initial stage of small carbon cluster formation, on 
the road towards the fullerene structure, when the 
sp 2 phase is energetically less favorable^ Furthermore, 
the interpretation of some IR bands from the inter- 
stellar dust could be related to the presence of sp- 
coordinated carbon. 5 Carbon chains are also interesting 
systems because they can be found connecting different 
graphene fragments. 6,7 Experimental work has been re- 
cently devoted to the production and investigation of 
sp-hybridized carbon structures4' 7 -'&^ i 10 ' 11 Furthermore, 
pure-carbon solid containing carbynoid structures have 
been produced by cluster beam deposition in the last 
years. 8 Carbynes survive landing because they are stabi- 
lized and protected by the sp 2 network. 

Open carbon (C) chains of any length must be ter- 
minated by molecular complexes to ensure stability, of 
course. In principle, linear carbon chains can be ob- 
tained with two different terminations: sp 3 termination, 
resulting in carbon atoms linked by alternated single 
and triple bonds (polyyne) and so with alternating bond 
length, and sp 2 termination, resulting in double bonds 
(cumulene) Experimentally polyynes are found to be 
more stable than cumulenes.— In this work, we have cho- 
sen hydrogen (H) terminations. The choice of saturating 
the chains with one or two hydrogen atoms reflects the ex- 
perimental fact, that the carbynes attached to graphene 
fragments can be bound with a variable number of carbon 
atoms. Optical spectra of carbyne are related to the typi- 
cal 7r-bonds of carbon valence electrons, which are degen- 
erate only in the polyyne case. Electronic spectra have 
been measured before for different polyyne molecules in 
neon matrices-^ or in gas phased 



The present article reports on the ab initio study of 
the optical properties of cumulene and polyyne proto- 
types, within the density-functional theor y 15 i 16 (DFT) 
using the independent-particle approximation— (IPA) in 
order to investigate if cumulenes and polyynes can be 
distinguished by means of their optical spectra. The 
calculations have been carried out using the ABINIT 18 
code. In order to support our findings, we have compared 
them with results within the Random-Phase approxima- 
tion (RPA) including local-field (LF) effects performed 
using the YAMBO code J*. We have derived our conclu- 
sions in two independent ways: on one hand side we have 
compared the low- and medium-energy spectral structure 
based on the DFT-IPA for our prototypes; on the other 
hand side we have drawn conclusion based on general 
spectral features which are obtained in both, the DFT- 
IPA and the RPA-LF results. 

Polyyne prototypes are simulated using short (seven or 
eight C) carbon chains saturating the ends with one H 
atom at each side, while cumulene prototypes are simu- 
lated with two H atoms terminating each side; the latter 
case is also studied for different orientations of the planes 
to which the terminating CH2 belong, and the differences 
in the spectra are reported. The terminating groups con- 
tribution has been singled out by applying the so-called 
real-space cutoff technique ; 20 ' 21 ^ 22 which was originally 
developed for surfaces and which has been adapted here 
for molecules. Using this technique, the influence of the 
H termination of the carbon chains on the optical spectra 
has been determined. 

This article is organized as follows: first, we intro- 
duce shortly the methodological background on which 
our calculations are based as well as the computational 
details (Sect. [IT]). Then we inspect the groundstate prop- 
erties of CyH^ and C§H„ chains (n — 2,4) with various 
lengths and H arrangements, where we investigate espe- 
cially their bonding character and their magnetization 
(Sect lTn)) . The bonding character and the spin polariza- 
tion are reflected in their electronic properties (Sect. IIV[) 
on which the final results, the absorption spectra, are 
based (Sect.|V|. In the latter section, we discuss the spec- 
tral structure obtained within DFT-IPA by performing a 
detailed peak analysis as well as inspecting more gen- 
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eral spectral features. Besides, the influence of LF to the 
DFT-IPA is discussed. Finally, we summarize and draw 
our conclusions. 



II. THEORETICAL BACKGROUND AND 
COMPUTATIONAL DETAILS 

In this section, the groundstate theory and the com- 
putational details are summarized as well as the basics 
for the calculation of optical spectra. 



A. Groundstate 

Our investigation starts from ab initio total energy cal- 
culations using the ABINI T 18 ! 23 code, which employs 
the density-functional theory scheme 15 (DFT) within 
the local-density approximatio n > 25 (LDA). The eigen- 
functions used to solve self consistently the Kohn-Sham 
equations 1 ^ are expanded in plane waves and, for the 
ionic cores, we have chosen norm-conserving pseudopo- 
tentials in the Troullier-Martins style.. 2 S ABINIT is de- 
signed for a periodic-cell approach. The molecule is 
placed in a supercell and surrounded by vacuum to pre- 
vent the coupling to its periodic images. Due to the 0- 
dimensionality of the system only the T point is necessary 
to sample the reciprocal space. Therefore, the conver- 
gence of the groundstate depends on the kinetic-energy 
cutoff which determines the number of plane waves in the 
expansion and on the size of the supercell. 

For an accurate groundstate, convergence requires a 
kinetic-energy cutoff of 20 Ha and a tetragonal supercell 
with lengths a = b = 25 ae and c = 40 ae- The car- 
bon chain is oriented along the z direction (c axis). This 
choice of the convergence parameters yields total energy 
differences less than I0~ 4 Ha. The carbynes have been re- 
laxed till the remaining forces are less than 5 TO -5 Ha/ 0,3. 
For details, see Ref. 1271 . 

YAMBOi^ works directly on the ABINIT ground- 
state results and no additional groundstate calculations 
are necessary. 

B. Optical properties 

The calculation of optical properties, in particular the 
photo-absorption cross section <j(lu) as a function of the 
energy oj is directly connected with the calculation of the 
polarizability function a, since^ 

a(uj) = — Im[47ra(w)] (1) 
co 

holds. Herein, Im[d(w)] is the imaginary part of the av- 
erage (over the x, y, and z) polarizability function and cq 
is the velocity of light. Thus it is sufficient to determine 
Im47ra(o;). 



In a first step, Im[47ra(u;)] has been calculated in the 
independent particle approach (IPA)jiS The probability 
Pvk,ck of the transitions between valence (v) and empty 
(c) states with electronic eigenenergies E v k and at a 
given k point in the reciprocal space can be calculated as 
the matrix elements of the velocity operator. With these 
probabilities, Im[47rav(a;)] (y — x,y,z) can be written 

Im[47ra„H] = ^^^^|p; kck | 2 

k v,c (,2) 

x5(AE-hu>) , 

where AE — E^ — E v ^, and m and e are the electronic 
mass and charge. In this approximation, only the matrix 
elements P^ k ck and the electronic eigenenergies are re- 
quired. For molecules, the sum over k can be reduced to 
the use of the T point only. Working within the DFT- 
IPA, local-field, self-energy, and excitonic effects are ne- 
glected. 

In order to determine the influence of hydrogen ter- 
mination of the chains, we have employed the real-space 
cutoff technique (see, e.g., Hogan et aL— ) to the DFT- 
IPA spectra. With this technique the contribution of a 
selected region of the simulation cell can be filtered out 
by including a boxcar function in the calculation of the 
matrix elements. This function is equal to one in the de- 
sired region and zero elsewhere. We have implemented 
the matrix elements with and without real-space cutoff 
to the ABINIT— code and we have tested the implemen- 
tation successfully.— 

Beside the DFT-IPA, we have employed more elabo- 
rated methods in the framework of the TDDFT in the 
linear response formalism^ In brief, this approach re- 
quires the solution of the Dyson equation for the den- 
sity response function. Here, also the standard Random- 
Phase approximation (RPA) with the inclusion of local- 
field (LF) effects is accessible. In order to obtain the po- 
larizability function for interacting particles (which is di- 
rectly related to the response function) one has to invert 
the Dyson equation. Usually, for periodic cell systems 
this is done in reciprocal space which can be quite cum- 
bersome for supercell systems. Therefore, an alternative 
method is provided by the Casida approach ) 33 i 34 i 35 i 36 i 37 
which moves the task to the eigenvalue problem 

with A = AE + U Q + f xc and B = U + / xc , where / xc 
is the exchange-correlation kernel and Uq is the Coulomb 
interaction term. For details see, e.g., Refs [36| and l32l. 
Eq. © is implemented in YAMBO±&2£ and allows the 
utilization of various approximations. Setting B = 
yields the Tamm-Dancoff approximation (TDA) whereas 
choosing / xc = the standard RPA-LF. 

We have performed calculations using YAMB 19 i 28 
employing the RPA-LF using both, the Dyson equation 
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FIG. 1: (color online) Sketch of the carbynes CgH2 (upper 
panel), CstU in the D2h (middle panel) and in the D^d con- 
figuration (lower panel). Note the different bonding type of 
the molecules. 



solution as well as the Casida approach. 



III. GROUNDSTATE PROPERTIES 

In this study we consider the carbynes CsEb and C8H4 
as representatives for even atomic polyynes and cumu- 
lenes (Cg chains). In the latter case, the terminating hy- 
drogens can have different orientations yielding the D 2 h 
and r>2d symmetry. The D 2 h has a planar configuration 
(all H and the carbon chain itself are in one plane) while 
for the T)ih one pair of H atoms is placed in a plane 
perpendicular to the one containing the other pair of H 
atoms, see Fig. Q] The same configurations are also stud- 
ied for chains containing seven carbon atoms (C7 chains) 
as odd atomic prototypes. Besides, also longer polyynes 
with an even number of carbon atoms are inspected. Here 
we will discuss the groundstate properties in terms of 
bond length and magnetization. 
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FIG. 2: (color online) Equilibrium bond lengths of C7 chains 
(upper panel) and Cs ones (lower panel). 



change whereas the average length of the single bonds is 
slightly reduced to values around 1.31 A. 

Inspecting C7H4 one obtains bond lengths very similar 
to the case of CsB^ due to the double-bond character of 
the C-C bindings. The situation is completely different 
for C7H2 due to the odd number of C atoms. In CsH2 the 
C1-C2 and C7-C8 bonds at the end of the chain are both 
triple bonds and the alternation of single and triple bonds 
continues for the whole chain due to the even number of C 
atoms. For C7H2, the C-C bonds at the end of the chain 
are triple bonds which is reflected in the bond length of 
1.23 A similar to the corresponding bond length of CsEb. 
The neighboring bonds have a length of 1.31 A which is 
between the length of the single and the double bonds 
of CsBi2. The chain "tries" to obtain alternating bond 
length. In the middle of the chain, the alternation from 
both ends gets in conflict and results in bonds with length 
of 1.26 A. 



B. Analysis of the magnetization 



A. Bond length 

The relaxation of C8H4 results in nearly equal C-C dis- 
tances of 1.27-1.28 A due to the double-bond character 
of all internal C-C bonds whereas for CsH2 the single and 
triple bonds result in bond length of 1.33 and 1.23 A, re- 
spectively (see Fig. [2]), as expected from the literature. 
The distances C1-C2 and C7-C8 are slightly different due 
to the presence of the terminating H atoms. As one can 
see in Fig. [2J the configuration {D 2 h or D 2 d) of CsH4 
does not influence the bond length significantly. Since 
the bond length of CsILi are equal, no change is expected 
for longer chains. Extending CsH2 to C12H2 and C16H2, 
one observes that the length of the triple bonds does not 



Analyzing the groundstate properties we obtained the 
expected result, that C8H4 in the ~D 2 d configuration is 
magnetic, whereas CsfLi D 2 h and CsH 2 are not. Since 
the T>2d differs only by the angle between the planes in 
which the H atoms are placed, we have investigated how 
the variation of angle between the planes containing the 
CH2 terminations affects the magnetization of the cu- 
mulenes, varying the angle from 0° (T>2h) to 90° (D2<j)- 
For this reason, we have performed total-energy calcula- 
tions forcing the Kohn-Sham orbitals to be occupied each 
with two electrons (spin up and spin down) like in the 
non-magnetic case of T) 2 h and on the other hand side, 
we imposed the spin configuration of the T) 2 d cumulene 
(magnetic). With these two setups, we have calculated 



4 




FIG. 3: (color online) Groundstate energy of CgH4 (upper 
panel) and C7H4 (lower panel) as a function of the relative 
angle between the H-H terminations. The solid line marks 
the magnetic (M) and the dashed line the non-magnetic (UM) 
configuration. 



the total energy as a function of the rotation angle. The 
result for the C 8 chains is presented in Fig. [3] (upper 
panel). Note that in the case of the D 2 ^ configuration 
the absence of spin polarization is energetically preferred 
by 59 mHa, whereas for T>2d the spin-polarized setup is 
more stable than the spin-unpolarized one by 11 mHa. 
Turning the H terminations, the system prefers to re- 
main in the unpolarized configuration until an angle of 
78°. After this angle the magnetic configuration becomes 
energetically advantageous, as we found for CsH 4 in the 
T)2d symmetry. These findings are in agreement with the 
results of Ravagnan et al£ performed with a different 
numerical implementation of the DFT equations. 

Interestingly, for the C7 chains the situation is opposite 
to the Cg ones. Here, the C7H2 and C7H4 with D 2 ^ sym- 
metry are magnetic, whereas the C7H4 in the D 2( j config- 
uration is non- magnetic. Performing an analysis similar 
to the C8H4 case, one finds that the D 2 h and the configu- 
rations until an angle of 5° are magnetic, while the config- 
urations for larger angles and the D 2( j are non-magnetic 
as shown in Fig. [3] (lower panel). Here, the energy differ- 
ence between the magnetic and the non-magnetic config- 
uration for T>2h is 4 mHa while it is 86 mHa for D 2t j. 



IV. ELECTRONIC PROPERTIES 

Optical transitions are excitations of electrons from a 
occupied state to an empty one in case of absorption. 
This is the basic assumption for DFT-IPA calculations. 
Furthermore, optical spectra can be interpreted inspect- 
ing the electronic states with respect to their (relative) 
energies and the charge density corresponding to these 
states. Thus, it is useful to inspect the electronic prop- 



erties of the carbynes investigated here in detail. 

The calculated electronic eigenenergies for our carbyne 
prototypes are schematically shown in Fig. [4] Since in 
the calculation of DFT-IPA spectra only energy differ- 
ences are considered, we have shifted the energy range in 
such a way that the highest occupied molecular orbital 
(HOMO) is at eV. In case of spin polarization we have 
chosen the HOMO of the major spin. The states close 
to the energy gap between the HOMO and the lowest 
unoccupied molecular orbital (LUMO) are energetically 
well distinguishable. Thus, one expects a clear structure 
of non-overlapping peaks in the low-energy range of the 
spectra. 

The electronic states of C 8 H 2 close to the gap are dou- 
ble degenerate due to the cylindrical symmetry of the 
chain. Because of this, the coupled ir x and ir y single 
orbitals are at the same energy. The triple and single 
bonds are clearly distinguishable in charge density of the 
HOMO state (see Fig.©. The linear sp bonds are present 
at lower energies, since they are more localized and feel 
more the attraction of the positive ions. LUMO states, 
instead, are placed around the odd bonds. There is a 
strong overlap between HOMO and LUMO around the 
carbons, in particular if we consider CI and C8. Because 
of this, one expects that a significant contribution to the 
optical spectrum can be attributed to electronic transi- 
tions in the spatial region of the last carbon atoms next 
to the H terminations. 

For C8H4 in the T>2h configuration the cylindrical sym- 
metry is broken, so the electronic states are no more de- 
generate, and there is an alternation of x and y polar- 
ized states. The HOMO of C 8 H 2 splits into C 8 H 4 (D 2 /,) 
HOMO-1 and LUMO, since the corresponding charge 
density is very similar. An analogous situation occur for 
the LUMO of C 8 H 2 (see Fig.©, which splits into HOMO 
and LUMO+1. The overlap between HOMO and LUMO 
is zero, since they are orthogonally polarized, so we do 
not expect an absorption peak for a transition involv- 
ing these states. On the contrary, HOMO and LUMO+1 
show a charge-density distribution on the same plane. 
The same applies to HOMO-1 and LUMO. For the lat- 
ter states, we see in Fig. [5] an important overlap around 
the termination, because these two states are polarized 
orthogonally to the hydrogen plane. The contribution of 
this region should be corrected by the real-space cutoff 
in order to eliminate the H contribution to the spectra. 

A mixed situation arises for C 8 H4 in the D 2£ ; symme- 
try. Here only the bands between HOMO-4 and LUMO 
are double degenerate for spin up and the bands between 
HOMO-3 and LUMO+1 for spin down. Majority and 
minority electron densities have a very similar shape, 
therefore only one component is shown in the Fig. [5] 
The 7r-derived HOMO-1 can be understood as a result 
of the merging of the HOMO-1 and HOMO of the planar 
conformer. 

Regarding the degeneracies and the splitting of eigen- 
states, the C 8 and the C7 chains are very similar. How- 
ever, the different bonding character of C7H 2 is reflected 
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FIG. 4: (color online) Schematic presentation of the rela- 
tive positions of the DFT eigenenergies for the considered 
C8 chains (upper panel) and the C7 ones (lower panel). For 
magnetic chains, the eigenenergies for the major (M) and mi- 
nor (m) spin is broken down. The energy scale is shifted 
to HOMO=0 eV (of the major spin for magnetic configu- 
rations). The main optical transitions within DFT-IPA are 
marked with dashed lines. 



in the charge-density distribution of the HOMO and the 
HOMO-1. While for CsH2 the charge density is placed 
around the triple bonds for the HOMO, for the C7H2 
it is localized at C atoms. Especially in the HOMO-1, 
the difference due to the even/odd number of C atoms 
can be seen: for C7H2 the corresponding charge density 
covering the three C atoms localized in the center of the 
molecule. 

Looking forward to the calculation of the optical spec- 
tra and the use of the real-space cutoff, the polyynes and 
the cumulenes should be treated differently. Cumulenes 
are chains of the style C„_2(CH 2 )2 and thus the termina- 
tion of the chain contains not only H but also C. Because 
of this, the cut should be placed between the last two C 
atoms of the chain and cutting of a complete CH 2 group. 
Inspecting the charge density for the states close to the 



FIG. 5: (color online) Charge density for Cs chains (left panel) 
and C7 ones (right panel); for each panel, from the fop to bot- 
tom, for the polyyne (upper panel), the T>2h cumulene (middle 
panel), and the D 2d one the HOMO-1, HOMO, LUMO, and 
LUMO+1 are presented from left to right. 



gap, only a little amount of charge is cut in the mid- 
dle and therefore, unphysical spectral features might be 
avoided. On the contrary, for the polyynes the termina- 
tion consists only of H resulting in a cut between the C 
and H atom. As shown in Fig. [5l for states close to the 
HOMO-LUMO gap, nearly no charge density would be 
cut. 

In general, the charge-density distribution is very regu- 
lar leading the conclusion that a similar distribution will 
be found even for longer chains. The only chain where 
charge-density differences might be found is C7H2 where 
the regimes showing alternating bonds and double bonds 
have different extents for longer chains. 



V. OPTICAL PROPERTIES 

In this section we present our results for the optical 
spectra of our polyyne and cumulene prototypes. In a 
first step, we analyze the spectral structure of our DFT- 
IPA spectra where we compare the peak positions in the 
low- and medium-energy range (Sect. IV A[) . Here, we 
have investigated also the contribution of the terminat- 
ing H atoms to the spectra using the real-space cutoff 
technique (Sect. |y"B) . 
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FIG. 6: (color online) Comparison of the imaginary part of 
the polarizability function as a function of the energy for the 
Cs chains (upper panel) and the C7 chains (lower panel). Sm 
correspond to major spin, S m to minor spin spectra. The 
peaks are labeled with A for polyynes, B for T>2d cumulenes, 
and C for D2/1 cumulenes. Here only the low-energy range is 
shown. 

Since the DFT-IPA yields fast but not highly accu- 
rate results, an alternative scheme for distinguishing the 
considered six prototypical chains has been developed 
based on general features observed in the optical spec- 
tra (Sect. |VC|) . In order to confirm the basic ingredients 
for this scheme, we have performed RPA-LF calculations 
for the non-magnetic systems C 8 H 2 , C 8 H 4 (D2/1), and 
C 7 H 4 (Dm) (Sect. ED). 



A. Spectra within DFT-IPA 

For the calculation of the imaginary part of the dielec- 
tric function it was possible to reduce the kinetic-energy 
cutoff to 18 Ha without loosing significant accuracy (for 
detail with respect to this parameter, see Ref. [3lh . This 
leads to a faster but equally accurate calculation. An- 
other important parameter is the number of empty states 
to be considered (see Eq. [5]) . For molecules the contin- 
uum states should be excluded, since an excitation of 
the molecule to these states corresponds to ionization. 
Continuum states are extended in space and thus, the 
eigenenergies corresponding to these states vary with the 
size of the supercell. An inspection of the energy of many 
states as a function of the size of the supercell yields the 
result that the energy of the states up to LUMO+7 is 
not influenced significantly by the size of the simulation 
box while the others are. Thus, we have considered in 
our calculation 8 empty states. Beside the parameters 
mentioned here, all other parameters are kept the same 
as for the groundstate calculations. 



The calculated spectra of the imaginary part of the po- 
larizability function as a function of the energy for C 8 H 2 , 
C 8 H 4 (D 2 h and D 2d ), C 7 H 2 , and C7H4 (D 2/l and D 2d ) 
are presented in Fig. [fj] (low-energy range) and Fig. 
(medium-energy range) , where also the resolution in ma- 
jor and minor spin contributions is given for magnetic 
compounds. 

Inspecting the low-energy range of the spectra of C 8 H 2 
(Fig. [6]), one notice immediately that the three considered 
prototypes are clearly distinguishable by their spectral 
structure: C 8 H 2 has only one strong peak (marked with 
Ai in the figure) at 2.90 eV which results from transi- 
tions between the HOMO and LUMO states. This tran- 
sition involves electrons in the tt x and n y molecular or- 
bitals (see Fig. [5j , which are degenerate due to cylindri- 
cal symmetry. Thus, the peak is polarized in z direction. 
For the non-magnetic C 8 H 4 D 2 /j cumulene the main peak 
is split into two peaks, whose energies are 2.07 eV and 
2.85 eV. Since the cylindrical symmetry is broken, both 
the HOMO and LUMO states of C 8 H 2 are split in two 
states for C 8 H 4 (see Fig. gj|. Here the HOMO-LUMO 
transition is forbidden because the respective wave func- 
tions lie in the x-z and y-z planes. The peaks are due 
to the transitions HOMO-l^LUMO (2.07 eV, Ci) and 
HOMO^LUMO+1 (2.85 eV, C 2 ). The C 8 H 4 D 2d shows 
also a double-peak structure, however, one peak is due 
to a major-spin (2.76 eV, Bmi) and the other due to a 
minor-spin (2.26 eV, B OT i) transition. Both transitions 
are HOMO-LUMO ones. From the low-energy range of 
the DFT-IPA spectra one can conclude that the C 8 chains 
can be distinguished in optical spectra, since the polyyne 
prototype shows only one main peak whereas the cumu- 
lenes show a double-peak structure. The two cumulenes 
can be distinguished in symmetry due to the fact that 
one is magnetic and the other not. 

Since the C 8 H 2 polyyne might have not exactly the 
cylindrical symmetry in the experiment, we have also cal- 
culated the groundstate and the spectrum of a non-linear 
configuration of C 8 H 2 where we altered the carbon po- 
sitions in order to create a zig-zag configuration shifting 
the carbons of the molecule (which is oriented along the 
z axis) in the x or y direction by a displacement of about 
the 5 % of bond lengths. In this way we have broken the 
cylindrical symmetry. Because of this, the main peak Ai 
is split in two close peaks at 2.85 eV and 2.80 eV. For 
peaks at higher energy we have found that the higher 
is the energy of the peak, the higher is its energy split- 
ting. Nevertheless, this split is much less than the ones 
of the cumulenes (Cj., C 2 and Bmi, B m i, respectively) 
and thus, the conclusion from above remains valid also 
in this case. 

Inspecting the spectra of the C7 chains, one observes 
a clear spectral structure as in the C 8 case. CyH 2 is 
magnetic and the HOMO-LUMO transition splits into 
two peaks at 3.30 eV (Ami, major spin) and at 2.71 eV 
(A m i, minor spin). Similarly, the non-magnetic C 7 H 4 
(D 2 d) shows only one peak at 2.74 eV (Bi) due to the 
HOMO-LUMO transition. This peak appeared split 
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for the magnetic CsH4 D 2f j cumulene. The magnetic 
C7H4 D 2 /j shows four transitions: for the major spin 
we have the HOMO-l^LUMO (2.74 eV, C M i) and the 
HOMO^LUMO+1 (3.62 eV, C M2 ) one and for the mi- 
nor spin a HOMO^LUMO (2.23 eV, C ml ) and a HOMO- 
l^LUMO+1 (2.99 eV, C TO i) one. As for the C 8 case, 
comparing the C7 chains one can conclude that the three 
prototypes can be distinguished since we have four peaks 
for C7H4 r>2h, two peaks for C7H2, and one peak for 
C7H4 (D 2(i ). 

Comparing the DFT-IPA spectra for the C7 and the 
Cg, chains, from the spectral structure all types of chains 
could be identified in experiments if their magnetization 
is known. The four-peak feature of C7H4 T>2h is unique. 
We have two magnetic double-peak features (CsITi D 2( j 
and C7H2) which are distinguishable due to the transi- 
tion energy. We have one unique non-magnetic double 
peak for CsELi D 2 /j. The only chains which cannot be 
distinguished unambiguously from the low-energy range 
of the DFT-IPA spectra are C 8 H 2 and C7H4 D 2d which 
are both non-magnetic and show only a single peak which 
differs in energy only about 0.1 eV. Note, that in the low- 
energy range all peaks shown in Fig. [5] are z polarized. 

Unfortunately, DFT-IPA spectra usually can give only 
a rough idea of what might be seen in the experiment. 
Comparing spectra performed with TDDFT methods or 
those which are based on solutions of the Bethe-Salpeter 
equation (for an overview of methods and results see, e.g., 
Ref. [29h . one observes not only a shift of the peaks but 
also a shift of the oscillator strength to higher energies 
which is not uniform. Nevertheless, the low-energy tran- 
sitions will contribute also to the corresponding peaks at 
TDDFT level, though, possibly at less extent. Regard- 
ing the shift of peak positions and oscillator strength to 
higher energies, the medium and high-energy range can 
be enhanced. Therefore, we inspect the medium-energy 
range of our DFT-IPA spectra, too. 

For the medium energy range up to 7 eV (see Fig. [7]) a 
clear spectral structure is absent as well as for the high- 
energy range (not shown here). In this range, the DFT- 
IPA spectra shows only a low intensity, which is of about 
0.5% — 2.5% of the intensity in the low-energy range. The 
spectral features involve transitions from states down to 
HOMO-5 to states up to LUMO+5. Inspecting the tran- 
sitions, no general trends have been found. Considering 
the polarization of the transitions, most of the peaks are 
z polarized. The remaining peaks have an x and/or y 
component only. Considering the magnetic compounds, 
one sees that some of the peaks for major and minor spin 
coincide in energy. For C8H4 D 2( j we have two of these 
peaks, B M2 / m2 and B m4 /B M 4- For C7H4 D 2/l it is only 
one, namely C M5 / m5 and for C 7 H 2 it is A M3/m4 . 



B. Spectra using the real-space cutoff 

The DFT-IPA technique gives us the possibility to de- 
termine the contribution of the terminating H atoms to 
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FIG. 7: (color online) Same as in Fig [6] for the medium-high 
energy range. 



the spectra by employing the real-space cutoff introduced 
by various groups, see e.g., Refs. I20II21II22I . As discussed 
in Sect. IIV1 the best choice of eliminating the H termi- 
nation would be to cut the H atoms only for polyynes 
and the terminating CH 2 group for cumulenes. However, 
the resulting chains would not be comparable in intensity 
due to the different number of remaining C atoms. Thus, 
we have decided to eliminate the terminating H atoms 
together with the neighboring C atom (resulting in the 
CH 2 group for cumulenes and a CH group for polyynes) 
by applying the real-space cutoff for the polyynes and 
cumulenes in the same range. 

Inspecting the low-energy range, we have found a low- 
ering of around 25% of the main peaks for the Cs chains 
applying the real-space cutoff. This can be explained by 
the fact that we have cut also one C atom at each end and 
thus, the contribution of the H atoms is very small. We 
can also conclude, that the intensity of the peaks is de- 
termined by the number of C atoms in the chain which is 
confirmed by comparing the intensities of the Cs and the 
C7 chains in Fig. [6l Thus, one expects for longer chains 
not only spectra containing more peaks (since there are 
more non-continuum states to consider) but also with a 
larger intensity. 



C. Conclusions from general considerations 

Often DFT-IPA spectra are not able to predict a quan- 
titative correct peak structure. Thus, we have analyzed 
our spectra also from a more general point of view in 
order to find an alternative way to draw conclusion. 
At higher energies mainly x and y polarized peaks are 
observed, whereas only few z polarized feature appear. 
Looking at the z polarized peaks only, we have found a 
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gap of at least 2 eV between the last z polarized peaks 
and the following high-energy z polarized feature. This 
inset of the gap is placed at around 8-10 eV and de- 
pends on the kind of the chain. Similarly, in the low- and 
medium-energy range up to about 9 eV only few peaks 
with x and/or y component are found. This brings us to 
the conclusion that the spectra can be divided in energy 
range: in the low and medium energy range the z com- 
ponent plays a major role whereas in the high-energy 
range the spectra is determined by x and y polarized 
peaks. This observation is useful for further considera- 
tions, since it allows the counting of z polarized peaks in 
the low-energy range till the gap inset. 

Inspecting other general features of the chains like sym- 
metry and magnetization, we are able to draw conclusion 
without detailed peak analysis. 

In a first step, the six chains can be divided into two 
groups, the magnetic chains CsH 4 (T>2d), C7H2, and 
C7H4 (D2/1) and the non-magnetic chains CsH 2 , C8H4 
(T>2h), and C7H4 (T)2d), which can be discriminated by 
the presence of a magnetic moment. 

The second step is guided by symmetry considera- 
tions: in each group there is only one chain which has 
no xy symmetry: the magnetic C8H4 (D 2 /j) and the non- 
magnetic C7H4 (D2h) where the x and y components of 
the spectra differ. For these chains one expects spectral 
changes (in the high-energy regime where x and y polar- 
ized features appear) for x and y polarized light. 

For the last step we exploit the fact that the spectra 
can be divided into two zones: the range up to 8-10 eV 
where the spectra shows mainly z polarization and the 
range beyond 8-10 eV, where the spectra has mainly x 
and y polarization. This range is separated by a gap of 
at least 2 eV from the next small z polarized peak. Since 
for the polyyne chains some states close to the HOMO- 
LUMO gap are degenerate, one expects less peaks for 
C 7 H 2 and C 8 H 2 than for C 8 H 4 (D 2d ) and C 7 H 4 (D 2d ) in 
the magnetic and non-magnetic case, respectively. Since 
the spectra in this range show mainly z polarized peaks, 
we are able to count these peaks till the inset of the gap: 
in the energy range up to about 10 eV we have three 
peaks for CsH 2 , 6 peaks for CgiLi (D 2 ^), 5 peaks for CgH4 
(D 2 d Sm), and 7 peaks for CsH 4 (D 2( j S m ). Similarly, 
there are 2 peaks for C7H 2 (Sm), 3 peaks for C7H 2 (S m ), 
4 peaks for C7H4 (D 2( i), 4 peaks for C7H4 (D 2 ^ Sm), 
and 6 peaks for C7H4 {T>2h S m ). For the comparison of 
C7H 2 with C8H4 (T*2d) (both magnetic) this means that 
we have 2/3 (Sm/S,„) against 5/7 (Sm/S™) or in total 
6 against 8, while for CsH 2 and C7H4 (D 2 ^) there are 3 
against 4 peaks. 

With these simple considerations we have shown, that 
the 6 chains can be distinguished by optical spectra with- 
out the necessity of a detailed peak analysis. 



VI. SPECTRA USING RPA-LF 

In order to inspect spectral changes using more ad- 
vanced methods, we have performed calculations using 
the RPA-LF. The main scope here is to confirm the va- 
lidity of the general conclusions, in particular, the exis- 
tence of the gap which is used to design a criterium for 
the peak counting and the number of peaks in the low- 
and medium energy range. For this reason, it was suffi- 
cient to do calculations not at full convergence, which is 
difficult to achieve. For this test we have performed cal- 
culations for the non-magnetic chains CsH 2 , CsH 4 (D 2 /j), 
and C7H4 (T>2d) only. 

The results show mainly the expected shift of peak 
positions and oscillator strength to higher energies. In 
detail, we have found a strong shift to higher energies of 
the peaks shown in Fig. [6l which now appear at energies 
around 5-6 eV, and an enhancement of the peaks with 
z polarization shown in Fig. [7] The shift is not uniform 
for the low-energy peaks. For example we have found 
for C8H4 only one strong peak followed by three small 
peaks whereas at DFT-IPA there are two strong peaks. 
Instead for C7H4 a double-peak structure with strong 
intensity appear followed by a less intense double-peak 
structure. It is possible that the first strong-intensity 
double feature merge to one peak at higher convergence. 
For CsH 2 a strong peak remained followed by two small 
ones. Nevertheless, the three non-magnetic chains still 
show a different spectral structure and are therefore dis- 
tinguishable by means of optical measurements. 

Regarding the general conclusions of Sect. IV CI we 
have found that for the z polarized peaks the oscillator 
strength are shifted from the high-energy peaks (above 
the gap) to the low-energy ones, while the x and y com- 
ponent remained nearly unchanged in intensity. Since the 
high-energy z polarized peaks have a very small intensity, 
a considerable spectral structure has been found only up 
to 8/9 eV followed by a gap of 2 eV. Thus, a counting 
of the peaks is also possible here. In particular, we have 
found 3 peaks for CsH 2 compared to 3 peaks DFT-IPA, 
four peaks for CsELi (D 2 /j) compared to 6 peaks DFT- 
IPA, and 5 peaks (2 strong double peaks followed by a 
tiny one) for C7H4 (D 2t j) compared to 4 peaks DFT-IPA. 
Also here, CsH 2 and C7H4 (D 2( z) are distinguishable by 
the number of z polarized peaks. 

In summary, employing the RPA-LF we find similar 
spectral features as in the DFT-IPA which have been 
used for our general conclusions. 

VII. SUMMARY AND CONCLUSION 

We have performed an ab initio investigation of car- 
byne prototypes inspecting their groundstate and optical 
properties in detail using DFT based methods. We have 
considered seven- and eight-atomic chains saturated with 
one H (polyyne prototype) and two H (cumulene proto- 
type) at each end, respectively. Polyynes and cumulenes 
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have a different bonding character. While the cumulenes 
show always double bonds, the polyynes have an alter- 
nating triple-single bond character. In case of odd-atomic 
chains, we have observed that this alternation results in 
the accidental building bonds similar to double bonds 
in the middle of the chain. Regarding the cumulenes, 
there are two possible symmetry configuration, namely 
T)2h and a D 2 d, corresponding to placing the saturating 
H atoms all in one plane or placing them in planes per- 
pendicular to each other. Depending on the symmetry, 
the C 8 H 4 D 2d is magnetic and C 8 H 4 D 2h is not. For C7H4 
the situation is reversed: here the T) 2 h chain is magnetic 
whereas C7H4 is non-magnetic. Rotating the H atoms, 
we have observed that the angle for which the cumulene 
becomes magnetic is different for even- and odd-atomic 
chains. Considering the optical spectra, we have found 
that the H terminations give only tiny contributions to 
the spectra. From the analysis of our DFT-IPA spectra 
we have seen that the six carbynes considered here are 
distinguishable by their optical spectra. However, due to 
well-known problems of the DFT-IPA a quantitative con- 
clusion cannot be drawn unambiguously. Nevertheless, 
we have been able to draw a scheme to distinguish the 
chains unambiguously based on the general features of 
the calculated spectra (within DFT-IPA and RPA-LF), 
the spin polarization of the systems, and their symme- 
try. With this scheme, the carbynes CsEb, C8H4 T) 2 h, 
CsH 4 D 2 d, C7H2, C7H4 D 2 h, and C7H4 D 2( z can be dis- 



tinguished. It might be possible to extend this scheme 
also for longer chains. 

To apply our scheme, one needs the knowledge of 
the magnetization of the system and polarization re- 
solved optical spectra are necessary. However, more ac- 
curate calculations might give quantitative spectra from 
which the chains can be distinguished by absorption 
spectra only without inspection of the magnetization of 
the chains and/or polarization resolution of the spectra. 
These calculations might be accelerated due to the fact, 
that the low-energy range of the spectra is mainly z po- 
larized and thus a detailed inspection of the x and y 
component is not necessary. We hope that this study 
here will encourage further theoretical and experimental 
investigations of carbynes. 
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